Attractive ultracold bosons in a necklace optical lattice 
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We study the ground state properties of the Bose-Hubbard model with attractive interactions on a 
M-site one-dimensional periodic — necklace-like — lattice, whose experimental realization in terms 
of ultracold atoms is promised by a recently proposed optical trapping scheme, as well as by the 
control over the atomic interactions and tunneling amplitudes granted by well-established optical 
techniques. We compare the properties of the quantum model to a semiclassical picture based on 
a number-conserving su(M) coherent state, which results into a set of modified discrete nonlinear 
Schrodinger equations. We show that, owing to the presence of a correction factor ensuing from 
number conservation, the ground-state solution to these equations provides a remarkably satisfactory 
description of its quantum counterpart not only — as expected — in the weak-interaction, superfluid 
regime, but even in the deeply quantum regime of large interactions and possibly small populations. 
In particular, we show that in this regime, the delocalized, Schrodinger-cat-like quantum ground 
state can be seen as a coherent quantum superposition of the localized, symmetry-breaking ground- 
state of the variational approach. We also show that, depending on the hopping to interaction ratio, 
three regimes can be recognized both in the semiclassical and quantum picture of the system. 

PACS numbers: 03.75.Lm 05.30.Jp, 03.65.Sq 31.15.Pf 



I. OVERVIEW 

Owing to the impressive progress in experimental tech- 
niques, ultracold neutral atoms trapped in optical lattices 
are nowadays widely recognized as a versatile toolbox 
bringing into reality ideal models of condensed matter 
physics [lj. A prominent example in this respect is no 
doubt the Bose-Hubbard model, originally introduced to 
sketch the physics of superfluid helium in porous media 
@ and subsequently shown to be realizable in terms of 
optically trapped ultracold Bosonic atoms 0, fl- This 
model, describing interacting Bosonic particles hopping 
across the M sites of a discrete structure, is characterized 
by a Hamiltonian of the form 
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A I 
m — l 



l (n m -l)-T ^ J > 



mm' Uj m u "m' 



(i) 



where a m , aL and n r 



are on-site Bosonic oper- 



ators, U measures the strength of the (on-site) boson- 
boson interaction, T is the hopping amplitude across 
neighboring sites and J is the so-called adjacency matrix, 
describing the lattice topology. Its generic entry, J mm ', 
equals 1 if the sites m and m' are adjacent, and oth- 
erwise. Parameters U and T are directly related to well 
defined experimental quantities, i.e. the scattering length 
of the Bosonic atoms and the intensity of the laser beams 
giving rise to the optical lattice, respectively 3]. The pos- 
sibility of tuning the lattice strength over a wide range of 
values played a fundamental role in the experimental ob- 
servation [^1 of the superfluid-insulator quantum phase 
transition predicted for Hamiltonian Q in the case of 
repulsive interaction pj. Further aspects of versatility 
of optically trapped ultracold atoms lie in the possibil- 



ity of tuning the atomic scattering length, and hence the 
boson-boson interaction, via Feshbach resonances H and 
realizing a great variety of different lattices, by means 
of a suitable laser beam setup 0, 0] , possibly involving 
holographic elements @ . 

In this paper we are concerned with the case of attrac- 
tive bosons in a ID periodic lattice, whose experimental 
realizability is made concrete by a recently proposed op- 
tical trapping scheme |£| and by the control over the sign 
and intensity of the boson-boson interaction granted by 
Feshbach resonances 0. This symmetric, necklace-like 
geometry is particularly interesting since, in the strong 
interaction limit, the system ground state is a quantum 
superposition of macroscopic spatially localized states, 
i.e. a so-called Schrodinger- cat state |9|- Although the 
state-of-art experimental setups ensure a remarkable con- 
trol and insulation from the environment, the realization 
of such deeply quantum states presents several difficul- 
ties, which will be briefly discussed in the final part of 
this work. Our main aim here is a study of the ground 
state properties of quantum Hamiltonian Q), which we 
compare to the corresponding state of a suitable semi- 
classical approximation to the Bose-Hubbard model. 

Before going into details, we mention that the attrac- 
tive Bose-Hubbard model has been addressed by two re- 
cent works ^3 El ■ In Ref. 01 a further nonlinear (cu- 
bic) repulsive term is introduced in Hamiltonian Q in 
order to prevent excessive localization for large values of 
the interaction. Such a term is ascribed to three body 
interactions coming into play at elevate atomic densities. 
However, experimental observations suggest that these 
interactions produce different phenomena resulting into 
a leak of condensate atoms rather than into a repulsive 
term [12|. Hence, here we do not consider such a re- 
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pulsive term, and assume that the total boson popula- 
tion is sufficiently low to rule out the above phenom- 
ena. In some respects our analysis is similar to that 
of Ref. whose results we compare to ours. Other 
recent results related to ours are found in Refs. [T^| . 
which also focus on annular geometries, yet in the ab- 
sence of a lattice potential. Finally, we mention that the 
attractive Bose-Hubbard model JIJ has received consid- 
erable attention as the quantum counterpart of the dis- 
crete self-trapping (DST) or quantum discrete non-linear 
Schrodinger (DNLS) equations 0, 0], particularly in 
relation to the issues of discrete quantum solitons [ltij 
and breathers [13, 0|- 

The layout of this paper is as follows. In Section [Q] 
we rapidly review the key features of the Bose-Hubbard 
model 0. Focusing on the special case of a periodic 
one-dimensional lattice comprising M sites, we recall 
the main properties arising from translational invari- 
ance. In Section II I II we describe a semiclassical varia- 
tional approach to Hamiltonian Q based on a number- 
conserving trial state. This approach results into a set of 
DNLS equations featuring a correction factor that makes 
them effective also in the deeply quantum regime char- 
acterized by large attractive interactions and small bo- 
son populations. Furthermore, we describe the three 
regimes that can be recognized for the ground state of the 
one-dimensional periodic lattice under investigation. A 
Schrodinger- cat regime for small values of the parameter 
t = T/\U\(Af— 1), where N is the total number of bosons 
in the system; a soliton regime, for intermediate values 
of t; a superfluid regime for large values of r. We show 
that the crossover between the Schrodinger-cat and soli- 
ton regimes occurs around a size-independent threshold 
T\ , while the critical point marking the boundary between 
the soliton and superfluid regimes, ? exhibits a clear 
dependence on the lattice size M, which becomes linear 
for large sizes. In Section HVI we illustrate the effective- 
ness of the semiclassical description introduced in Sec- 
tion [ffi] In particular, we show that in the Schrodingc- 
cat regime, the delocalized ground state of the quantum 
model is very well described by a coherent superposition 
of the degenerate localized ground states of the semiclas- 
sical picture. Furthermore, we discuss the relevance of 
the semiclassical thresholds t% and t^ 1 ^ in the quantum 
picture. In particular we show that the local observables 
considered in Ref. are not sensitive to size dependent 

threshold , which is related to the finiteness of the 
system, but rather to the size-independent threshold ti, 
which is related to the discrete nature of the system. We 
also show that T2 is signaled by nonlocal observables. 
In Appendix ^ we show that the variational trial state 
introduced in Section II I II is actually a su(Af) coherent 
state, which is a further justification (?) for the effective- 
ness of the semiclassical picture. In Appendix iBl we show 
that in the Schrodinger-cat regime the localized ground 
state of the semiclassical picture coincides — at the first 
perturbative order in r — with a coherent superposition 



of the almost degenerate M lowest states of Hamiltonian 

id- 
ii. HOMOGENEOUS PERIODIC ID LATTICE 

Let us now summarize the main properties of Hamilto- 
nian ||TJ. First of all, we recall that the relevant Hilbcrt 
space is in principle infinite, and it is spanned by the Fock 
vectors generated by the action of the positive powers of 
the creation operators aj^ onto the vacuum state |fi), i.e. 
the state such that a m |f2) = for any m = 1,2, . . . , M . 
However, Hamiltonian commutes with the total num- 
ber operator, N = Y] m —i n m , and hence its study can be 
restricted to the finite size eigenspaces of N. Henceforth, 
we assume to be in the eigenspace relevant to the total 
population N, whose size is d(Af,M) = As 
we mention, we are interested in the case of attractive 
bosons (U < 0) on a ID homogeneous periodic lattice. 
Hence we particularize the adjacency matrix in Eq. 
to 

where, on a periodic lattice comprising M sites, the site 
labels and M + 1 correspond to sites M and 1, respec- 
tively. With this choice Hamiltonian is clearly trans- 
lationally invariant. Hence all of the relevant observables 
cannot break this symmetry, i.e. exhibit a dependence 
on the site label. Owing to the translational invariance 
the eigenstates of H can be labeled by a further quan- 
tum number. This is made quite evident by switching to 
the reciprocal space, i.e. by introducing the transformed 
operators 

M 

b k = M- 1 ' 2 e l2 * krn/M a m , (3) 

ra— 1 

making the kinetic term of Hamiltonian (JTJ diagonal. Of 
course, the reciprocal space has the same size d(Af, M) 
as the direct space, and it is spanned by the reciprocal 
Fock vectors, generated by the action of the powers of 
the bp's on Note indeed that N = Yl P =i bpbp and 
that b p \Q,) = for any p. In this representation the dis- 
placement operator changing a m into a m+ i is simply 

2tt M 
£ = exp(z— P), P = Y J Pb\b P . (4) 

P =i 

Since [D, H] — 0, the Hamiltonian can be studied sep- 
arately in the M eigenspaces of D. Note indeed, that 
D M = I, which entails that the eigenvalues of the dis- 
placement operator are the M th roots of unity. After 
noticing that the reciprocal Fock vectors are clearly eigen- 
states of D, the block diagonal structure of H can be 
made evident by simply arranging these states accord- 
ing to the relevant eigenvalue, which we denote by T>. 
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A perhaps more convenient choice for the displacement 
quantum number is the integer k € [0, M — 1] such that 
— ilog(X>) = 2nk/M. In particular, the ground state of 
the system is the lowest eigenstate of the block relevant 
to T> = 1 or, cquivalently, k = 0. Introducing a further 
integer label j ranking the eigenstates of of Hamiltonian 
(|T|) according to their energy, we set 



D\1> 



0)\ 



with < iy+i- In particular, \ipi'') and E\"' denote 
the ground state of H and its energy. We recall that this 
state is translationally invariant, while, more in general, 
all of the eigenstates of H are delocalized. That is to 
say, these eigenvalues cannot give rise to site-dependent 
expectation values of local operators, such as e.g. n m . 

This means that, at least in principle, the states \ip[ k>> } 
exhibit the features of the celebrated Schrddinger cat 
state. Indeed, as it is well known, for T <C U the ground 
states of the M blocks of H are almost degenerate coher- 
ent superpositions of macroscopically localized states of 
the form 



e M 



0)\ 



WW 



(5) 
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r (0) 



M 
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■\Af m ), Wm) = 



t \M 



(<) 

VAT! 



■in). (6) 



macroscopic local phase and population, identified with 
arg(</> m ) and \<fi m \ 2 ■ The request that the trial state 
satisfies the Schrodinger equation on the average 
(Q\idt — H\&) = 0, leads to the identification of S with 
the effective Lagrangian for the dynamical variables 4> mi 
and 



H 



M-1U ^ 



Af 2 ^ 

m—1 



\<f> m \ 4 -T 



M 

E 

m.m' — l 



J, 



mm' tym' 



(8) 



with the corresponding effective Hamiltonian |23|. Min- 
imizing TL with respect to the </i m 's provides a semiclas- 
sical variational approximation to the ground state of 
the Bose-Hubbard Hamiltonian. Introducing the chem- 
ical potential fi, i.e. a Lagrange multiplier allowing for 
the normalization constraint in Eq. Q , the ensuing equa- 
tions are 



Af-l 
Af 



U\(j) m \ 2 4>m — 



mm' H>m' 



(9) 



We recall that Eqs. 10 are actually the fixed point equa- 
tions corresponding to the equations of motion for the 



This can be quite easily seen by treating the kinetic term 
of Hamiltonian JQ) perturbatively [l^. In Sec. IIVI we 
show that the Schrddinger cat structure of Eq. JHJ can 
be extended to a whole interval of T/U ratios, by replac- 
ing the Fock vectors \Af m ) with suitable su(M) coherent 
states which are strictly related to the localized ground 
state of the semiclassical approximation to Hamiltonian 



III. TIME DEPENDENT VARIATIONAL 
PRINCIPLE 



As we recall above, the attractive Bose-Hubbard 
model Q has received significant attention also as 
the quantum counterpart [l4( of the so-called discrete 
self-trapping (DST) or discrete non-linear Schrodinger 
(DNLS) equations [lflL l2?j . l2ll | . These semiclassical equa- 



tions are equivalently obtained by discretizing the Gross- 
Pitaevskii equations j2^| or applying to Hamiltonian 
a time dependent variational principle (TDVP) based on 
Glauber's coherent states [2^. Here we adopt the same 
approach as Ref. [2j|, using a macroscopic trial state of 
the form |$) = e tS \$) where 



i 



M 



v j v .j v \m — 1 



in), E l^l 2 =^ ( 7 ) 



m—1 



instead of a product of Glauber's coherent states. The 
complex quantity 4>m appearing in Eq. Q is a dynamical 
variable describing the bosons at lattice site m through a 



^E J « 



(10) 



obtained via the canonical Poisson brackets {<j) mi <j)* m ,} = 
ifim m' ■ Here by fixed points we mean the "normal modes" 
whose time evolution is entirely contained in a simple 
phase factor e~ lfJ,t . Equations l|ll)fl clearly belong to the 
DNLS class and, as we show in the following section, they 
are able to capture many of the features of the quan- 
tum ground state also in the extreme quantum regime 
of large interactions and small populations. This is due 
to the extra (JV — 1)/M factor appearing in their in- 
teraction term, which is absent in the DNLS equations 
standardly adopted in the semiclassical approach to the 
Bose-Hubbard model. The presence of this extra factor 
is a direct consequence of the macroscopic trial state, 
Eq. J7J, being an eigenstate of the total number opera- 
tor N . This is not true of the Glauber coherent state 
product adopted as a trial function in standard TDVP, 
which conserves the total number of bosons only on the 
average j23. Hence the resulting dynamical equations, 
which are perfectly equivalent to standard DNLS equa- 
tions, do not feature the extra (Af — 1)/A/" factor. The 
effectiveness of the description provided by the macro- 
scopic state (7J| can be explained by noticing that it is 
actually a coherent state itself, of su(M). Note that Af 
is the algebra representation index of the present two- 
boson realization of su(M) as confirmed by the fact that 
N commutes with the algebra generators a^a m / for any 
1 < m,m' < M. Hence the conservation of the boson 
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population is already incorporated in this coherent state 
picture . This is illustrated in Appendix lAl 

As it is well known, the nonlinear terms can break the 
translational symmetry characterizing Hamiltonian (JHJ) 
with the choice in Eq. Indeed, if the interaction U 
is sufficiently strong, the ground state of the system is 
arbitrarily localized at one of the M equivalent lattice 
sites. This ground state is a fixed point of Eqs. iJTUjl. and 
can be easily determined numerically by treating Eqs. @ 
as a set of self-consistent linear equations (2^. More in 
general, depending on the value of the T/U ratio, it is 
possible to recognize three regimes for the semiclassical 
ground state. For large ratios, the ground state is super- 
fluid, the boson population being evenly spread over the 
entire lattice. For sufficiently low values of the T/U ra- 
tio this delocalized supcrfluid solution becomes unstable, 
and the ground state of the system becomes soliton-likc, 
featuring a localization peak breaking the translational 
symmetry of Hamiltonian (jSJ . The critical point marking 
the boundary between the superfluid and soliton regime 



IS T 



AM) 



[2Msin 2 (7r/M)l 



we define r = T/\U\(Af - 1) HH IM |2j 



)] * where we recall that 

[H 0,113. As t is fur- 
ther decreased, the localization peak — which initially 
involves a significant number of lattice sites — becomes 
sharper and sharper, until the entire boson population 
is basically concentrated at a single lattice site. This 
situation is the is the semiclassical, symmetry-breaking 
counterpart of the Schrodinger-cat state, where the bo- 
son localization is symmetrically distributed over all of 
the lattice sites. As we discuss in the sequel, it is pos- 
sible to recognize a threshold, or rather a crossover, be- 
tween this Schrodinger-cat regime and the soliton regime 
where the localization peak involves a large portion of 
the lattice. Note indeed that in the latter situation 
the periodic boundary conditions and the discretization 
of the lattice produce a small effect, and the ground 
state is well described by the solution to the nonlinear 
Schrodingcr equations on a ID infinite domain, <f> m — 
e^ v /7V/2Asech(TO/A), where A = 4T/\U\(J\f- 1) = 4r 
is the characteristic width of the population peak. The 
inset of Fig. ^ shows that this is actually a satisfac- 
tory approximation. Hence, the crossover between the 
Schrodinger cat and soliton regime can be estimated by 
requiring that this characteristic width equals the lattice 
spacing, i.e. Awl, which results into r = t\ « 1/4. 
Note that this is equivalently obtained by equalling the 
energies of the system in the two extreme situations 
T <C \U\ and T > \U\. Indeed in the former case the 
ground state is completely localized and its energy is 
Af (Af — l)U/2, while in the latter case the ground state is 
delocalized and its energy is —2TAf. We emphasize that 
exactly the same estimate is obtained considering the en- 
ergy of the quantum ground state. Conversely, the der- 
ealization transition at r = tv^ ^ , marking the boundary 
between the soliton and superfluid regime, corresponds 
to the peak width being of the same order of the lat- 
tice size, A M. Note indeed that ri M ^ is linear in M 




FIG. 1: (Color online) Normalized interaction energy, 
Eq. for two different lattice sizes, M — 20 (light solid 

line) and M — 40 (dark solid line). The derealization tran- 
sition clearly occurs at r = t% ■ Furthermore, the plotted 
quantity exhibits a size independent behavior in the small r's 
regime. The crossover at r = n corresponds to a change in 
the convexity plotted function. This is related to the ultravi- 
olet cutoff due to the discreteness of the lattice. Note indeed 
that in the continuum limit the corresponding quantity would 
diverge as r -1 for r — ► 0. The inset shows a comparison be- 
tween the localized ground state of Ti (bullets) on a lattice of 
size M — 40 and the corresponding solution in the continuum 
limit (solid line). 



above are quite evidently signalled by the behavior of the 
(normalized) interaction energy 



E 



M 



Af 2 



M 
Af 2 



ill) 



for sufficiently large lattices 



The regimes discussed 



as a function of r, where (•) denotes expectation value on 
the trial state J7J equipped with the ground state solution 
to the fixed point equations J^J- This is shown in Fig.^ 
where two lattice size are considered, M — 20 and M = 
40. Recalling that (n m ) — Af/M, it is easy to see that 
the quantity in Eq. is also related to the fluctuation 
of the local density of bosons. 



IV. SEMICLASSICAL VS. QUANTUM 
PICTURE 

In this section we compare the features of the semiclas- 
sical picture introduced in the preceding Section to the 
corresponding features of the quantum model, Eq. l|T]l. 
First, we direct our attention to the energy and structure 
of the ground state. Then, we discuss the relevance of the 
semiclassical thresholds in the quantum picture. We con- 
sider lattice comprising up to M = 20 sites, and popula- 
tions consisting of up to TV = 10 bosons. All of the quan- 
tum observables are obtained evaluating the eigenstates 
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FIG. 2: (Color online) Relative energy difference between 
the exact quantum ground state and the relevant semiclas- 
sical approximation. The curves refer to boson population 
M = 10 and lattice sizes M = 6, 8, 10, 12. The inset shows 
the energy range spanned by the ground states E[ k ^ (shaded 
region) and the first excited states (region between the 
dotted lines) of the blocks relevant to different values of the 
quantum number k associated to the translational invariance 
of H [see Eq. ©]. These result refer to M = 12 and N = 10. 
The vertical dashed lines signal the crossover at t = t\. 



of Hamiltonian Q and the relevant energies by means 
of Lanczos-like algorithms. The translational symmetry 
of Hamiltonian Q is conveniently taken into account in 
order to reduce significantly the computational demand. 



A. Ground-state energy and structure 

As it can be seen in Fig. [21 the ground state energy of 
the quantum model, Eq. is in general remarkably well 
described by the relevant semiclassical result. Indeed the 
maximum in the relative error settles at about 2% as the 
lattice size increases. Furthermore, the agreement be- 
tween the quantum and the classical result is extremely 
satisfactory both in the large r and in the small r region. 
The agreement in the former region is expected, since in 
the noninteracting limit the quantum ground state of the 
system is exactly given by Eq. J7J) with <p m = y/M/M, 
which is nothing but the superfluid semiclassical ground 
state taking over after the derealization transition at 
T ( M ) _ The agreement in the strong interaction limit is 
more surprising, and it requires further discussion. In- 
deed, in this region the semiclassical and the quantum 
ground states are significantly different, the former being 
strongly localized at one lattice site, while the latter being 
always completely delocalized. As we recall in Sec.|H] the 
eigenvalues of the quantum Hamiltonian Q are labeled 
by a quantum number k ensuing from translational in- 
variance, Eq. JSJ). The inset of Fig. |3 clearly shows that, 




FIG. 3: (Color online) Normalized boson density ($j\n m \$j) 
as provided by the quantum localized state (solid lines), com- 
pared to the corresponding semiclassical quantity \4>m\ 2 (bul- 
lets) for N = 10 bosons on M — 12 lattice sites. 



in the small r region the actual ground state, ), is 
almost degenerate with the ground states of all of the 
remaining Hamiltonian blocks, Hence, in this re- 

gion, the actual ground state of the system is satisfacto- 
rily described by any superposition of the states \ip[ k ^). 
Of course this is not an eigenstate of H , but the relevant 
dephasing time is extremely long, being inversely propor- 
tional to the energy spreading of the involved eigenstates. 
In particular, one can choose a superposition of the form 

H3 

M 

l^-HEe-^SfclV^) (12) 
fc=i 

which features a population peak localized at lattice site 
j. We emphasize that the phase factors jk are not arbi- 
trary, but have to be chosen to make sure that the \tp[ k>> ys 
are phase coherent. Recall indeed that a generic eigen- 
state is defined up to an arbitrary phase factor. The 
superposition in Eq. i|12[) features the desired localiza- 
tion properties only if the 7fc's compensate for possible 
phase differences among the {ip^^'s 29]. As it is shown 
in Fig. |21 the boson density (<f>j\n m \$j) is in remarkable 
agreement with the corresponding semiclassical quantity, 
\4> m \ 2 for sufficiently small values of r. At larger values 
the agreement is lost since the energy spreading of the 
eigenstates involved in Eq. I|12l) becomes significant, and 
the superposition cannot be confused with the ground 
state any more. Furthermore, as it is shown in the inset 
of Fig. |2 the energies of the \ip^Ys become comparable 

to those of the \ij>\ 's. 

By resorting to perturbation theory, it is possible to 
show that the agreement between the trial state (Q) and 
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the localized state Eq. (fT^|) . goes beyond the lo- 
cal density of bosons. Indeed, at the first order in the 
perturbative parameter r, they are actually the same 
state. This is discussed in Appendix [5] Notice that 
the above discussion entails that, in the strong coupling, 
Schroedinger-cat regime, the actual ground state of the 
quantum Hamiltonian Q is a coherent superposition of 
localized states of the form {7J), where the </> m 's are ob- 
tained solving the corresponding semiclassical problem 
ifffil. Once again, we emphasize that the extra (A/ - — 1)/W 
in the DNLS equations (jlOf) . stemming from the choice 
of a number-conserving trial state, Eq. J7J), plays a fun- 
damental role in ensuring the remarkable agreement be- 
tween the quantum and the semiclassical picture, Eqs. 
Q and JBJl, respectively. 



B. Quantum signature of the semiclassical regimes 



Let us now discuss the relevance of the semiclassical 
thresholds t\ and rj*^ in the quantum picture of the 
system. This issue has been considered to some extent 
in Ref. ^l|, where only the derealization transition oc- 
curring at r = Trj, is taken into account. The authors 
report that this critical point is signalled by the behavior 
of some quantum observables, such as the one-particle 
and A/"-particle correlation functions, 



61 = —r(a 



l n a m+l) 7 



M , f . M M 



(13) 



where (•) denotes expectation value over the ground 
state. More precisely, they observe that, with increasing 
interaction, C\ features a sudden drop around the semi- 
classical derealization threshold r% , while CV features 
a surge roughly at the same point. Furthermore they re- 
port that this threshold is also signalled by the number 
distribution, i.e. the probability p v that the boson pop- 
ulation at a given lattice site equals v. More precisely, 
they report that this distribution becomes double peaked 

t ^ ( M ) 
for t < 7-3 . 

Here we show that the above quantum observables do 
not signal the derealization threshold, but rather the 
crossover between the soliton and Schrodinger-cat regime 
occurring around r = t\. First of all, we observe that 
the variation of these observables is not so sudden to 
identify a precise transition point. They rather signal a 
crossover interval. Furthermore, in the case considered in 
Ref. [ll|, M = 6, N = 10, the two semiclassical thresh- 
olds n = 1/4 and = 1/3 are relatively close, and 
both compatible with such crossover interval. However, 
as we illustrate in Section ITTT1 t\ = 1/4 is substantially 

size-independent, while t^ 1 ^ clearly depends on the lat- 
tice size. A study of the scaling properties of the above 
quantum observables shows that the features considered 
in Ref. do not depend on the lattice size, and occur 
in the region around n |30l |. 




I 



o M=12,A^=10 

+ M=6,A'=10 

M=20,A'=5 

* M=10,A'=5 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
X 

FIG. 4: (Color online) Top: contour plot of the number 
density p v {r) for Af = 10 bosons on a lattice of size M = 12. 
The darker the shade, the lower the density. The white bullets 
signal the value »„(t) at which the density is minimum. The 
vertical lines signal the semiclassical thresholds. Bottom: plot 
of Vui(t) for several choices of M and M. 



Let us first of all consider the number distribution p v , 
i.e. the quantity such that (n m ) = Y^Lo^Pf = A/"/M 
and, in general, (ri^) = Y^f—p v 9 Pv It is clear that in 
the strong interaction limit this distribution features two 
peaks at v = and v = JV, while it is substantially 
vanishing for all of the remaining numbers |llj| . Recall 
indeed that in this limit, , i.e. deep in the Schrodinger-cat 
regime, the ground state has the form in Eq. @ llLll9|. 
In this regime the number fluctuations are very large. 
Actually, as we discuss in Sec. 11111 this situation is the 
quantum counterpart of the semiclassical situation where 
the boson population is almost entirely concentrated at 
one of the lattice sites. Hence, the ground state cannot 
be influenced by the boundary conditions or by the lat- 
tice size, and this regime is expectedly entered at r as t\. 
This is actually confirmed by Fig. which illustrates 
some key features of the number distribution. As it can 
be seen in the contour plot of in the top panel, for small 
r the number distribution p v (r) features two maxima at 
the extreme values v = and v = M, while its minimum 
is at v — Af/2. As soon as t becomes larger than t\ the 
relative maximum at v = M disappears, and the mini- 
mum suddenly moves to this extreme value. The location 
of such minimum, is signalled by the white bullets on 
the contour plot. Note that no special feature can be rec- 
ognized in correspondence of the derealization threshold, 
signalled by the vertical dashed line. The bottom panel 
of Fig. 0] reports the (normalized) location v m {r)/N of 
the minimum of Pu{t) for several choices of the lattice 
size M and of the boson population N . In all cases the 
sudden displacement of the minimum signals that the 
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ealization thresholds relevant to the lattice sites considered 
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number distribution becomes single peaked at the semi- 
classical crossover ri, i.e. at the crossover between the 
Schrodinger-cat and soliton regime. Basically the same 
discussion as above applies to the A/'-particle correlation 
function CV defined in Eq. (|13fl . ft is indeed clear that it 
must exhibit a significant surge when the Schrodinger-cat 
regime is entered and the number distribution becomes 
double peaked, i.e. around r = t\. Indeed, for this corre- 
lation to be non-vanishing it is necessary that the prob- 
ability of finding all of the N bosons at the same lattice 
site is significant. This is confirmed by Fig. reporting 
the behavior of CV(t) and of its derivative with respect 
to t for several choices of the boson population and of 
the lattice size. No particular feature occurs at the M 
dependent derealization thresholds, signalled by the ver- 
tical dashed lines. Actually, a dependence on the lattice 
size M is hardly recognizable, the curves referring to the 
same number of bosons Af being almost perfectly coinci- 
dent. Conversely, the location where the plotted quantity 
becomes significant and its concavity changes are satis- 
factorily estimated by the semiclassical size-independent 
crossover at r = t\. Figure [|J] shows that also the sub- 
stantial drop in the single-particle correlation function 
C\ considered in Ref. [ll| actually occurs around r = t\ . 
This is confirmed by the behavior of the derivative of this 
quantity, shown in the bottom panel of the same figure. 
Again, the plotted quantity is not sensitive to the size- 
dependent delocalization threshold Tq , essentially due 
to the fact that it is a local quantity, as it can be seen 
from its definition, Eq. (|13fl . 

The size-dependent delocalization threshold between 
the soliton and superfluid regime is signalled by a typ- 
ically nonlocal quantity such as the occupation number 



b bo of the lowest reciprocal mode, where bk is defined in 
Eq. J3J). This is clearly evident in Fig. [7J where (&J& )/A/" 

is plotted with respect to both r (top panel) and t/t^ 1 ^ 
(bottom panel), for several lattice sizes and boson popu- 
lations. Unlike the quantities in Figs. I4lb1 {b\bQ)/N' ex- 
hibits a clear dependence on the lattice size, when plotted 
against r. Conversely, when the abscissae of the plot are 
rescaled so that all of the size-dependent delocalization 
thresholds occur at the same point (signalled by the verti- 
cal dashed line), a very satisfactory data collapse becomes 
evident. This means that the drop in the occupation of 
the lowest reciprocal mode is related to the delocalization 
transition observed in the semiclassical model. 



V. OUTLOOK 

In this paper we consider the Bose-Hubbard Hamilto- 
nian describing ultracold attractive bosons trapped 
into the necklace-like lattice whose realization in terms 
of optical potentials is described in Ref. ||. In particu- 
lar, we study the properties of the ground-state of such 
extremely symmetric quantum Hamiltonian, which in the 
strong interaction limit is known to exhibit the features 
of a Schrodinger-cat state. We introduce a variational 
approximation to this ground state based on a su(Af) 
coherent state, Eq. {7J). The equations governing the 
dynamics of the labels of this coherent state belong to 
the class of DNLS equations, yet they feature an unusual 
correction term ensuing from the total number conserva- 
tion. Owing to this correction term, the ground-state of 
this semiclassical approximation provides a remarkably 
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satisfactory description of its quantum counterpart not 
only — as expected — in the small-interaction superfluid 
regime, but also in the deeply quantum Schroedinger-cat 
regime, i.e. for strong attractive interactions. In particu- 
lar we show that in this regime the delocalized quantum 
ground state of the system can be seen as a coherent 
superposition of the localized ground states of the vari- 
ational approach. Furthermore, we show that the three 
regimes clearly recognizable in the semiclassical picture 
are signalled by suitably chosen quantum observables. 
The latter result is compared with similar findings re- 
ported in Ref. [Tl| . In the intermediate regime the de- 
scription provided by the above semiclassical picture is 
not as satisfactory as in the two outermost regimes. The 
structure of the quantum ground state in this parame- 
ter range is different, and probably much more complex, 
than that of the variational trial state Q, and requires 
further investigations. 

Before concluding this paper we briefly consider some 
experimental issues. First of all, we observe that attain- 
ing the deeply quantum Schrodinger-cat regime might 
prove very challenging from the experimental point of 
view. Indeed, this would be most likely obtained by start- 
ing from the ground state of the (almost) noninteracting 
superfluid regime, which is relatively easy to realize. By 
making use of Feshbach resonances and/or by increasing 
the strength of the optical lattice, the ratio T/\U\ can 
be subsequently decreased down to suitably small val- 
ues, so that the system is driven into the Schrodinger-cat 
regime. Note that, independent of the value of such ra- 
tio, the ground state of the system has always the same 



symmetry, i.e. it always belongs to the same eigenspace 
of the displacement operator defined in Eq. Q. This 
means that, at least in principle, only the relevant block 
of Hamiltonian is involved in the process. Actually, 
this is true provided that the system is perfectly isolated 
from the environment. Under this rather ideal assump- 
tion, the achievement of the Schrodinger-cat ground state 
basically depends on the adiabaticity of the parameter 
variation. However, the interaction with the environment 
much likely introduces a coupling between the ideally un- 
coupled blocks of Hamiltonian Q, which has dramatic 
effects in the strong interaction regime. Indeed in this sit- 
uation the actual ground state is almost degenerate with 
the lowest states of the blocks relevant to k = 2, . . . , M, 
and an arbitrarily small perturbation can induce a tran- 
sition to one of these excited states. 

Another issue worth mentioning is the actual possibil- 
ity of making measurements on a system such as the one 
considered here. As we recall in Section^ there is an up- 
per bound to the number of attractive bosons M that can 
be loaded into the optical lattice. Indeed, the maximum 
attainable boson density (n m ) = Af must be sufficiently 
low to rule out the complex phenomena described in Ref. 
|12|. On a one-dimensional system like the one consid- 
ered here, the maximum boson population compatible 
with such constraint can be too small to allow measure- 
ments, , such as e.g. absorption imaging. This difficulty 
can be possibly circumvented by considering a stack of 
coaxial necklace optical lattices isolated from each other, 
in place of a single necklace || . Assuming that these lat- 
tices can be loaded with the same number of bosons, they 
are expected to behave like copies of the system discussed 
in this paper. On the other hand, the larger overall boson 
population should allow for effective measurements. 
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APPENDIX A: SU(A/) COHERENT STATES 

In this section we give a sketch of the proof that the 
variational trial state introduced in Eq. J7J) is actually a 
su(M) coherent state, where M is the number of sites 
composing the lattice. According to the standard theory 
[3ll 13^ , the coherent states of a given group G are rep- 
resented, up to a phase factor, as |cs) = <?|ext), where 
g G G and |ext) is a so-called extremal vector. The latter 
is identified by determining the maximal isotropy sub al- 
gebra B € Q c , where Q c is the complexified version of 
the algebra Q associated to the group G. 

In our particular case, Q is the su(M) algebra rep- 
resented in terms of the M lattice boson operators 
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a m . Hence Q — {h s ,a\ n a m i : m ^ m'; m,m' = 
1,2, ... , M; s = 1, 2, . . . , M — 1} and a generic element 
of the SU(M) group is given by 
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where the Cartan operators h s are suitable linear com- 
binations of the on-site number operators n m = a^am, 
x m G K and cr mm i G C. The relevant maximal isotropy 
sub algebra B is formed by (complex valued) linear com- 
binations of all of the generators of G, except the M — 1 
operators of the form a^dm, where m corresponds to one 
of the lattice sites, and m ^ fh. The extremal state |ext), 
is defined as the state such that 6|ext) = Ab|ext) for any 
b G Q n B, with A b G C 32]. This yields |ext) = \J\f & ), 
where \Afm) is defined in Eq. ©. Of course, all of the M 
possible choices for m, and hence for the extremal state, 
are equivalent. For definitcncss, here we set m = 1. 

It is possible to show that for any <? G G there exists 
an unique decomposition g = g z e b , where b £ Q n B and 
the so-called displacement operator has the form 



exp 



M 

Ef 

m=2 



* t 



(Al) 



Hence, up to an irrelevant phase factor arising from e b , 
one gets |cs) = g|ext) = g z \Afi) 

On the other hand, the variational trial state intro- 
duced in Eq. Q can be rewritten as 



(^) 



M 



i M 

|0), A=-=^2<p m a. 



(A2) 



rn — 1 



Observing that A f = Ba\B\ and hence that (A^)^ 
B{a\) N B\ where 



EM - M n / T T \ 

^jj, — e cos (7, ^ — le Uk , u — l^m=2 u m> 



one gets 



|<J> =B^Lr\Sl)=B\Af 1 ), 



(A3) 



where we made use of the fact that B^\fl) = |f2). After 
some further algebraic manipulations, state l|A3|) reduces 
to (up to an irrelevant phase factor) 

|$\ _ e iyiA/' e Em=2(CmaJ„Ol+COlOm) \J\f \ ^ 

where £ m = 8 m e lipm G C. The comparison of the lat- 
ter expression with Eq. (|A1|) fully discloses the SU(M)- 
coherent-state structure of the variational trial state |<J>) 
defined in Eq. Q. 

It is important to notice that, while |$) defined as in 
Eq. I|A2(I is operationally useful to derive the semiclas- 
sical picture relevant to Eqs. JSJ) and I|1U|) . the hidden, 



equivalent form of |<I>) as a SU(M) coherent-state (CS) 
indeed is responsible for such a semiclassical scenario. 

In general, CS are known to provide a formulation of 
quantum- mechanical problems [3l|, |33] particularly close 
to the corresponding classical Hamiltonian picture. Let 
us consider a classical Hamiltonian He for 2M degrees 
of freedom and the corresponding quantum model Hq 
obtained via canonical quantization. The dynamics of 
the quantum model can be investigated enacting a time- 
dependent variational principle (TDVP) based on a co- 
herent state |S), with 3 = (£i,...,£m) G C m . This 
allows one to recast the Schrodinger equation for Hq 
into a classical problem whose Hamiltonian is 7i(S) = 
(S|.Hq|S). The corresponding Hamiltonian equations de- 
scribe the evolution of the CS labels £ m , which thus play 
the role of dynamical parameters. Simple analytic for- 
mulas show that £ m 's can be rewritten in terms of the 
canonical variables appearing in He (see, e. g., Ref. 
[23} ~). More in general, the fact that |S) is a coherent 
state ensures that the domain where 2 ranges (CS mani- 
fold) coincides with the phase space where the dynamics 
of He is represented. 

The CS variational method exhibits two further notice- 
able features. First, the classical Poisson brackets char- 
acterizing the S dynamics preserve the algebraic struc- 
ture of the operators constituting Hq and generating the 
group G. Also, the coherent-state labels £ m , beyond their 
canonical role in the TL dynamics, are the expectation val- 
ues of Hamiltonian operators that embody information 
on microscopic interaction processes. For example, in the 
specific case under investigation, ($|a+ afc|$) = 
($|n m |<fr) = 1 4> m 1 2 in Hamiltonian JSJ arc manifestly re- 
lated to a+afe and n m in Hamiltonian (Q. Notice that 
Hamiltonian (jHJ matches exactly the classical version of 
Eq. m as soon as N / (N — 1) ~ 1. An interesting 
results in the present work is the profound similarity be- 
tween the (purely quantum) properties exhibited by the 
ground state of Hamiltonian (JIJ and the properties of 
its semiclassical minimum-energy configuration derived 
from state l|A2() in the Schrodinger-cat regime. 



APPENDIX B: PERTURB ATIVE COMPARISON 
OF THE QUANTUM AND SEMICLASSICAL 
GROUND STATES 



In this section we show that in the strong interaction 
limit the superposition defined in Eq. I|12|) coincides with 
the coherent state where the m 's correspond to the 
semiclassical solution to Eqs. © localized at lattice site 
j. We do this resorting to perturbative theory, i.e. by 
assuming that 



M 

E 



i) 



(Bl) 
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and 

M 

V = -T Jmm>al n a m ' (B2) 

m,m' — 1 

are the unperturbed and perturbative terms of the quan- 
tum Hamiltonian (JIJ, respectively. At the th order in 
t = T/U(N — 1) the quantum superposition |<Ej) in 
Eq. [fl"2"|l coincides with the Fock state \J\fj) defined in 
Eq. ©. The latter belongs to the M-fold degenerate 
lowest eigenspace of Hq. Indeed, 

HalAfj) = E \Afj), E Q = ^N{N - 1) (B3) 

As to the first-order perturbative correction to it is 
standardly obtained as |<E,)i = (H — E )~ 1 V\Af 3 ) . After 
some simple algebraic manipulations we get 

|*i>i=rV^(|^+) + |^)) (B4) 

where 

\r,+ ) = ±Jl == ^\si) Wit) = ^=^=L\{1). 



Hence, to the first order in r, 

m*m+rVtf(\rit) + \ri7)) (B5) 

On the other hand it is easy to show that, to the first 
perturbative order in r, the solution to the semiclassical 
equations 10 localized at lattice site j is given by 

<j) m ~ VAf 8 m j + tVN (Smj+l + iSmj-l). (B6) 
Plugging this expression into Eq. (7J gives 

m^l^+rVMdvp + lrj-)), (B7) 

i.e. the same result as in Eq. I|B5(1 . Hence, to the first 
perturbative order, the quantum superposition of the M 
lowest eigenstates of H defined in Eq. I|12|) coincides with 
the state obtained plugging the suitable localized solution 
to the semiclassical equations (JHJl into the coherent state 
0. We emphasize that in the strong interaction limit 
the zero temperature properties of H can be investigated 
based on the ground state of the corresponding semiclas- 
sical model, Eq. JSJ, which is evaluated at a considerably 
lower computational cost. 
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